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^ ■ 1 Abstract 

O 

• The recent success of a proposal from some time ago to explain the spectrum 

of cosmic rays allows some strong conclusions to be made on the physics of 
O ' supernovae: In the context of this specific proposal to explain the origin of 

-4— > ' cosmic rays, the mechanism for exploding supernovae of high mass has to be 

, the one proposed by Bisnovatyi-Kogan more than 30 years ago, which was then 

based on a broader suggestion by Kardashev: A combination of the effects of 
rotation and magnetic fields explodes the star. Interestingly, this step then 
leads inevitably to some further suggestions, useful perhaps for the study of 



gamma ray bursts and the search of a bright standard candle in cosmology. 

2 Introduction 

Cosmic Ray physics has inspired us for many years, and has given rise to many 
interesting books and reviews, some of which are [101 [T71 071 ITS] . 

Several burning questions exist in high energy astrophysics: 



What is the physics of the supernova explosion of massive stars? The most 
common idea to be explored over the last few decades is that the burst of 
neutrinos, certainly sufficient in energy, as shown by the supernova 1987A, 
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is doing it. However, the detailed mechanism has yet to be worked out 
successfully. 

• What is the origin of Galactic cosmic rays? What makes it possible for 
cosmic rays from Galactic sources to reach energies such as 3 10^*^ eV? 

• What is the mechanism of Gamma Ray Bursts, and what is it's relation 
to supernovae? 

• Is there a brighter standard candle possibly available for cosmology than 
the now so famous Supernovae type la? 

Here, in this review, we wish to combine several recent advances made to 
give a tentative answer to all these questions, starting with the work on cosmic 
rays: 



3 The cosmic ray spectrum 

We wish to explain the entire spectrum of cosmic rays, and focus here on Galac- 
tic cosmic rays, so the energies up to about 3 10^ GeV. Some time ago we made 
a proposal to explain Galactic and extragalactic cosmic rays in six consecutive 
papers, [101 1^ IHSl EZl EH] • The basic idea was to distinguish the different 
sites of supernovae into the lower mass stars, that explode into the interstellar 
medium, and the higher mass stars that explode into their own stellar wind. 
The basic data to explain are as follows, |96l l68]: 

• The overall spectrum is about approx. until the "knee" , which is a 
bend downwards at around 3 10^^ eV 

• approx. E""^'^ beyond the knee 

• a slight downward dip from 3 10^ eV, sometimes referred to as the "second 
knee" 

• a transition near 3 10^^ eV, with then approx. E~'^''^ again 

• uncertainty beyond 5 10^^ eV, either a mild cutoff (HIRES) or a continu- 
ation (AGASA) 

and in more detail, considering electrons and nuclei: 

• electrons E'"^-^ as well (from radio, low E), up to about 20 GeV 

• electrons E~^'^ observed to 3 TeV: 
loss dominated, so injection E~'^'^ 

• positron fraction a few percent 

• abundances enriched 

• anti-proton fraction about 10~^ 

• heavy elements tend to have a slightly flatter spectrum 



4 The arguments 

Already in 1934 the original suggestion by W. Baade & F. Zwicky was that the 
most energetic particles probably come from the most energetic phenomonea 
known to us at the time, supernovae, ^Hl- E. Fermi then provided an argument 
on how reflection at magnetic irregularities can enhance the energy of a charged 
particle, 1949 and 1954, |34[ I35j . The detailed physics of the Fermi-acceleration 
was worked out in a simple approximation in a series of papers, starting in 1977 
by I. Axford et ai, G. Krymsky, j^, R. Bell in two papers, ^Bl, and then again 
by R. Blandford &: J. Ostriker, |HJ with a nice extensive review by L. Drury, 



4.1 Supernova explosions into the ISM 

Then P. Lagage & C. Cesarsky, inO], worked out the details for the maximum 
energy of a particle subject to Fermi acceleration in the shock of a supernova, 
exploding into the interstellar medium. Presciently, they already used the self- 
similar expansion into a very tenuous medium, now known to exist from X-ray 
observations. For the details of the acceleration of energetic particles at the 
shock they use the concept that most of the acceleration happens in a regime, 
where the magnetic field is parallel to the shock normal - as well its disturbances 
- and provides most of the scattering. Here the expansion solution is a similarity 
or Sedov-solution, a self-similarity, where the outer shock radius r behaves as 
r ~ t^/^, and the shock speed is r ~ ^--3/2^ slows down considerably 

over time. This slow late evolution makes these supernova remnants dominant 
in observations, they just live a rather long time. 

• Emax ^ 10^ ZS_6Gey, using already the hot ISM of density 3 10-3 cm-3, 
confirmed only in 1997, |84j . The expansion is self-similar, usually referred 
to as a Sedov expansion, and yields then a maximum particle energy of 
about Z 10^'^ eV, an order of magnitude below the knee. 

• Maybe secondary acceleration is then required in the ISM - maybe by 
Shockwaves inside superbubbles, the motley mixture of powerful stellar 
winds, and young supernova explosions. This picture has mostly been 
explored by I. Axford. A quantitative prediction is still outstanding for 
this scenario. 



4.2 The knee 

The key problem, recognized immediately, is that this maximum energy does 
not even come to the energy at the knee, and so that feature cannot be explained 
in such an approach. Therefore many attempts have been made: 
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4.2.1 Transport? 

Cosmic rays get injected with a spectrum, which is subsequently modified by 
losses from the disk, and so obviously, to get a kink in the spectrum may mean 
a change in behaviour of these diffusive losses; such a kink could be visible 
as the "knee". This is a question of the transport of cosmic rays through a 
turbulent magnetic ionized gas; general discussions of the magnetic field are in, 
e.g., inHEHI- Such a change in transport implies a special scale in the interstellar 
medium, and would also lead to an increased anisotropy at the kink. However, 
there is no such scale as yet recognized. 

4.2.2 Turbulence? 

The evidence from the interstellar medium and from plasma physics suggests - if 
described by an isotropic approximation with an eye towards transport of cosmic 
rays - that the spectrum is a Kolmogorov spectrum. There are three arguments 
in favour of using a Kolmogorov spectrum for the context of cosmic rays: First, 
the difference in behaviour between nuclei, and electrons, allows the exponent 
of the diffusive law to be limited by about 0.4, so the energy dependence of 
the diffusion coefficient k is weaker than k ~ E^''^, Second, any steeper 

dependence would produce anisotropics at the higher energies, which are not 
seen, I^U]. And, finally, any steeper dependence would seriously dilute the 
number of stars contributing at the higher energies, making any reasonable 
powerlaw fit to the data meaningless, despite the fact, that such fits work well 
to describe the data, ^511^ . 

4.2.3 Spallation? 

Cosmic rays interact, and in fact most isotopes of Lithium, Beryllium, and Bor, 
as well as the odd-Z elements, and the sub-iron elements arise from cosmic ray 
interaction, the destruction of nuclei in a spallation process, see M. Garcia- 
Munoz et al., |381 139j . Therefore, a change in this interaction may lead to a 
loss of nuclei at higher energy; the interstellar medium is discussed, with some 
emphasis on the Galactic Center region, e.g., in JSB]. However, when working 
through this idea, one finds that the column for interaction must be much 
higher than what is observed; in fact, from the abundances of all the spallation 
products, as well as positrons and anti-protons we have a fairly good idea what 
the interaction actually is, and it is not enough to provide such a drastic change. 
Furthermore, that data suggest today, that the "knee" is a feature in rigidity 
E/Z, and so plasma physical processes are implied. 

4.2.4 New source? 

Obviously, in some proposals we require a new source at 3 10^^ eV, or maybe 
even at lO^"^ eV, and so one might realistically ask, whether a new source might 
not be necessary already at lower energies, if that source is supposed to work so 
well at higher energies. A new source is easiest to accomodate, if the new source. 



dominant at higher energy, has a flatter spectrum than the source dominant 
at lower energy; in such a case the transition is trivial, since the higher energy 
source just overtakes the lower energy source naturally in flux. However, if the 
transition is supposed to occur at or near the knee, then the higher energy source 
must have a steeper spectrum rather than a flatter spectrum, and so we require 
a multitude of strange coincidences: The higher energy source must have a 
peaked spectrum, peaking just when the lower energy source peters out, and so 
have the same flux at that energy, while not contributing at lower energy in any 
significant way. That is really hard to believe. So, it is then actually simpler, 
and has been suggested, that we use a new source from the start, and, e.g., 
adopt the point of view, that Gamma Ray Bursts produce all the cosmic rays, 
in work by C. Dermer (discussed in a talk at the Aspen meeting, January 2002). 
However, then we are back where we started, we still need an argument why we 
have a change in slope at a specific rigidity, the knee. Also, the expected flux 
of cosmic rays from Gamma Ray Bursts is nowhere near to what is necessary, 
[7^1 [THj . However, Gamma Ray Bursts should produce cosmic rays, and the 
question is whether there is any way to ascertain their contribution, perhaps 
through high energy neutrinos. Another possible source could be the activity 
in the Galactic Center, [611 164j . However, there is no observational evidence at 
this time, that such a localized activity could be dominant at the location of 
the Solar system. 

4.2.5 Change in acceleration? 

Therefore, what remains, is more mundane perhaps, a small modification in 
acceleration of the existing cosmic rays. The proposal by I. Axford et al. 
is that normal supernovae provide cosmic rays to about the knee, and that 
then in the supernova wind bubbles around the really massive stars, the cosmic 
rays get further acceleration, and then have a slightly diminished efficiency in 
acceleration; after all, the detailed physical properties are somewhat different in 
the environment of many supernovae and their predecessor winds, as compared 
to the interstellar medium. The only step missing in this very nice argument 
is a quantitative prediction as to what expect in this picture; one may surmise 
that the predictions may not turn out to be quite so different from the proposal 
described below. 

5 Proposal for the Origin of Galactic Cosmic Rays 

The riddle, which provoked a new approach now some time ago, was the key 
observation, that the radio shells of supernova remnants show a radial direction 
for the magnetic field, parallel to the shock-normal. This is in exact contra- 
diction to the simple plasma physics expectation that the enhancement of the 
magnetic field component parallel to the shock surface should make that com- 
ponent dominant behind the shock, see j22j . 

Additional support for a new approach came from the observation that radio 
features could be seen in motion, with an amplitude of velocity quite close to 
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the overall expansion speed of the supernova shock front, and yet measured 
with respect to the expanding frame of reference. 

Third, the thickness of the observed shells both in the radio and also to some 
degrees in X-rays, is much too large to allow a simple explanation. There must 
be some reason for these shells to be almost 1/3 of the outer radius oberved. 

Finally, there appeared a theoretical argument, that cosmic ray dominated 
shocks are unstable, in work by R. Ratkiewicz, I. Axford, G. Webb, e.g., iTQ^. 

Therefore, we introduced the key premise that the particle transport in the 
shock region is dominated by large scale turbulence; in such a case there is no 
time for any turbulent cascade, and so no time for a Kolmogorov cascade either. 
Therefore the transport is just given by a specific scale, and a specific velocity, 
naturally the thickness of the shocked region, and the velocity difference across 
the shock. This is in the spirit of the early work on turbulence such as by 
L. Prandtl and Th.v. Karman. This premise allows then to explain all these 
observed features. 

This means that the shock region is strongly turbulent, and the front of the 
shock, as seen along the line of sight at the edge of some observed supernova 
remnannt, is actually a superposition of many smaller shocks, advancing and 
retreating, like the waves on a beach. The geometric scale of these shocks 
is again given by the radial overall shock thickness, just given by momentum 
conservation. Therefore, curvature drifts do play a role, and introduce a scale 
in energy per charge EjZ. 

It can be surmised, that the limit cycle arguments by M. Malkov & L. O'C. 
Drury will eventually lead to a similar conclusion, j62j . Such a convergence of 
theoretical concepts will be an important test of the picture. 

5.1 Supernovae into a stellar wind 

The concept then leads to the consideration of supernovae exploding into their 
own stellar wind, j9,S | I2()| [^. Stellar winds may be considered in a first approx- 
imation as a Parker wind, with ~ sin 0/r in spherical coordinates. The dom- 
inant magnetic field character is an Archimedian spiral, with the field mostly 
very nearly tangential. The key aspects are i) that the wind is magnetic, with 
the magnetic field given by the star, and not the interstellar medium, and ii) 
that the wind has a density gradient ~ r~^, which means that the shock stays 
fast, and does not slow down as in an explosion into a medium of constant 
density. Also, magnetic fields in massive stars exist, they have been detected 
through nonthermal radio emission and maser emission, |21j . 

As the magnetic fields get very weak and also radial towards the pole of 
a Parker-type wind, one needs to consider shock acceleration in two limiting 
regimes: 

At the pole and near to it we need to use the acceleration dominated by the 
magnetic field parallel to the shock normal, as done by ^60^. However, at the 
equator, where the shock direction is perpendicular to the dominant magnetic 
field, we need to use another limit, as shown by R. Jokipii, jSl]. Therefore we 
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also need a matching condition, and this leads to a proposal for identifying the 
knee as arising with this matching energy, the maximum energy near the pole 
and close to it: 

Emax = ZerB{2,/A Vsh/cf (1) 
The maximum energy in the equatorial region, and in fact, most of 47r is 

E^ax = ZerB (2) 

Since B is inversely proportional to radius r in the wind, these energies are 
actually constant throughout the wind, see e.g. [H2I- We still need to show, 
that the particles actually have enough time to reach these maximum energies, 
but this appears plausible. We also need to show, why this leads to a different 
spectrum beyond the knee: This is due to the effect that curvature drifts, which 
provide a small proportion of the particle's energy gain, when subject to many 
shock transitions, get slightly weaker when the Larmor radius of the particle 
moves past the scale corresponding to the maximum energy (really E/Z) near 
the pole, see [231 for a very detailed discussion. 

Those stars which have powerful winds, are usually enriched in chemical 
elements such as Helium, Carbon and Oxygen, and so this line of reasoning 
actually picks up a theme, that Wolf Rayet stars are a major contributor to 
cosmic rays, |83| . Therefore, supernova explosions into stellar winds are pro- 
posed to provide most of the Galactic cosmic rays, with the notable exception 
of the element Hydrogen, for which they provide only some fraction. 



6 Stellar evolution 

For the physics of cosmic rays and supernovae we need to consider stellar evolu- 
tion in some detail. What happens to stars of different zero-age main sequence 
masses M: 

• Single stars with M < 8 Mq give no supernovae. 

• Stars in the mass range 8 < M < 15Mq explode as supernovae, but the 
explosion goes into the interstellar medium. Such explosions lead a Sedov 
type expansion, which becomes slow and has a long-lasting remnant. 

• Stars with 15 < M < 30Mq explode as supernovae with a substan- 
tial wind, but enriched only in Helium. The mass in the wind-shell is 
moderate. We refer to these stars as Red Super Giant stars (RSG). 

• Stars > 30 Mq explode as supernovae into a strong wind; the chemical 
abundances are strongly enriched, since the mass loss eats deeply down 
back into the star. The mass in the wind-shell is large. We refer to these 
stars as Wolf- Rayet stars (WR). 

Obviously, these numbers for the masses are approximations. 
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7 Predictions 

The predictions, all made some time ago, here focussing on wind-SNe, were 
described in [101 HH HH HI IM| • 

They gave the following quantitative spectra: 

• £'-2.74±o.04 £qj. ^j^g interstellar medium supernovae, with the abundances 
that correspond to the medium through which the supernova shock races. 
This spectrum is predicted to have a cutoff in the 100 TeV range. 

• ^-2.67-o.02±o.02 |-,giQ-^ knee for the wind supernovae, for which the shock 
races through the predecessor wind, and correspondingly has enriched 
abundances. 

• ^~3.07-o.07±o.07 above knee for those same stars winds. 

• A bending of the wind spectra at the "knee" near Z 10^^ eV, and a final 
cutoff at the "ankle" near Z 3 10^"^ eV. These specific rigidities have an 
uncertainty of about a factor of 2. 

• When the supernova explodes, a powerful shock wave races through the 
wind, and then smashes into the shell; this wind-shell is made up of both 
wind-material and interstellar medium material, from the environment of 
the predecessor star. Such a shell has two shocks at its boundaries, one 
slowing down the wind, and one speeding up the interstellar material. If 
there is no convective instability, then a contact discontinuity separates 
wind material from environmental interstellar medium material. 

• The shock is loaded with energetic particles from cosmic ray acceleration, 
and then these cosmic rays suffer from spallation in the shell, [271 128j . 
Both cosmic ray particles as well as wind material particles get broken 
up. This, in the case of WR shells, gives the Be, B, Li nuclei in cosmic 
rays, which serve then as tracers of cosmic ray lifetimes, and interaction. 
It also gives the odd-Z elements, as well as the sub-Fe elements in cosmic 
rays. 

• As a new prediction following from the work in [2H^ we find a small spectral 
fiattening is expected due to differential spallation for Fe-like nuclei. The 
fraction of nuclei removed by spallation from the set of Fe-like nuclei is so 
large, that across the range of energy per particle measured, this fraction 
is larger for lower energy, and so the spectrum flattened. This proposal 
is quantitatively consistent with the data. 

At that stage we did not differentiate the cosmic ray acceleration scenario 
for red supergiant winds (RSG), and Wolf Rayet (WR) winds. We now do 
differentiate those two kinds of stars with winds in their interaction in the 
wind shell, |2H]- The RSG stars have a weaker wind, therefore a lower mass 
shell as compared to the WR stars, and so we argued that the cosmic ray 
interaction, when the supernova shock finally hits the shell, is convective for 



RSG stars, and diffusive for WR stars. In this picture the RSG stars produce 
the gamma ray emission from the Galactic plane, and the WR stars produce the 
spallation products such as Beryllium, and give only some small contribution 
to the Galactic gamma ray emission at moderate photon energies. 

7.1 Implications 

Such a model implies, as had been predicted already by B. Peters in the fifties, 
[711 172j . that the turnover at the knee, and at the final cutoff is gradual, as 
the elements roll off in sequence of Z. Furthermore, since the knee and final 
cutoff energy are given by stellar parameters, the magnetic field, implicitly the 
rotation (the wind must be asymptotic already at the surface, which implies 
slow rotation at the surface) and the explosion energy, these parameters should 
be very nearly the same for all stars. 

This corresponds exactly to what has been argued in a very different con- 
text already by G. Bisnovatyi-Kogan in (1970), who based his argument in 
turn on a broader suggestion by N. Kardashev (1964), j^HEni- In this picture, 
the predecessor star is rotating and has a magnetic field. When the core of 
the star collapses, the core is spun up to a fast rotating disk, and the collapse 
along the plane of symmetry stops. The potential energy at that stage then is 
transmitted to the envelope by the torque of the magnetic fields (i.e. their angu- 
lar momentum transport), and that potential energy constitutes the explosion 
energy. 

Therefore, if the data are really well fit by the model, the conclusion is 
strongly suggested that not only the concept of magneto-rotational explosion 
works, but that in fact it works with very similar numbers in all stars. 

8 New data 

Since the original detailed quantitative predictions were made, new data have 
appeared, which provide very serious tests, e.g. |5T] : 

• At present the data set from the KASCADE experiment is the best across 
the knee region, and provides a first challenging test. As we show below 
they are very well fit by the data. Other data also show consistency, 
but with large error bars, and with unknown systematics still, so further 
conclusions can not be drawn from them. 

• In paper |^ we predicted the energy dependence of the B/C ratio, from 
interaction in the thick WR-star shells, as E~^f'^ . A refined treatment was 
given in [28| . As V. Ptuskin then showed some time later, at the ICRC in 
Salt Lake City, |74j . the best fit to the data gives an energy dependence 
of which is quite consistent. 

• In 1997 St. Hunter et al., [^, and M. Mori, |f)7j . showed that the gamma 
ray emission from the Galactic plane could only be fitted by a cosmic 
ray spectrum, which was substantially flatter than the one observed. In 



10 



P71 05] we suggested that the spectrum could be explained as interaction 
in the source region of cosmic rays, specifically the RSG star shells. The 
interaction happens then with a spectrum of E^"^/^ . We will discuss the 
details of this fit elsewhere. 

8.1 Specific tests 

The tasks ahead of us are the following, to be quite specific: 

• We normalize all abundances at TeV, using the data collection of |96j . 

• We need to remember, that Hydrogen has two components, only one of 
which may have a knee (wind-SNe, RSG), the other one has a cutoff below 
the knee (ISM-SNe); on the other hand Hydrogen is extremely low in WR 
winds. We need to note here, that the data do not force us at this time 
to assume that the knee of the RSG stars and the WR stars is the same; 
in fact, it could be, that both knee and cutoff for RSG stars are quite 
different, much lower in rigidity than for WR stars. In our first data fits, 
described below, we do not yet distinguish RSG and WR stars as regards 
the knee and final cutoff. The KASCADE data are forceful only for the 
abundances from the WR stars. 

• Also Helium may have two different observable knees or cutoffs, corre- 
sponding again to RSG supergiants, and WR supergiants, but they also 
might be same; we have to keep this uncertainty in mind. 

• We then extrapolate with the predicted spectrum to the knee in rigidity, 
bend at the knee in rigidity, and extrapolate to the cutoff. As the cutoff is 
also in rigidity, there is a rolloff of the overall spectrum, starting from the 
cutoff in Helium, which we can actually identify with the "second knee". 

• We then compare to KASCADE, Akeno, HIRES, ... One important point 
is to remember, that any fit to the data across the knee has to "reach" 
the data at higher energy, such as from Akeno, etc. This restricts any 
possible fits quite severely, see (5T] . 

8.2 Work to be done 

The work is spread among many of our partners, as follows: 

• The detailed abundances are being considered with A. Popescu (Bukarest), 
and N. Langer (Utrecht). 

• The energy in explosions is considered with G. Pavalas (Bukarest). 

• The positron spectra are done with W. Rhode et a/.(Wuppertal). 

• The gamma ray spectrum of the Galaxy is done with S. Casanova (Bonn), 
R. Engel (Karlsruhe), W. Rhode (Wuppertal), and many others. 
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• The anti-proton spectra are done with E.-S. Seo, R. Sina (Univ. Mary- 
land), and R. Engel (Karlsruhe). 

• The detailed fit to the KASCADE data is first described in this paper, in 
the following main section. 

• Later we will also consider a more detailed fit to the Akeno, AGASA and 
HIRES near the ankle with S. Ter-Antonyan (coauthor here), A. Vasile 
(coauthor here), Zh. Cao (Utah) and St. Westerhoff (New York). 

9 Cosmic ray properties from air shower data 

The following section was written by S. Ter-Antonyan. 

In general, the relation between energy spectra (SQ^/S-Eq) of primary nu- 
clei (A) and measured EAS size spectra at observation level {AI / AN*^^) is 
determined by an integral equation 



where the kernel function (dW/dN*^^) of the equation depends on primary en- 
ergy, type of primary nucleus {A = 1,4, .. . ,59), zenith angle (9), A — Amt 
interaction model and response functions of measurements. 
Eq. © is a typical ill-posed problem and has an infinite set of solutions for 
unknown primary energy spectra. However, the integral Eq. © turns to a 
Fredholm equation, if the type of a primary nucleus is defined directly in the 
experiment (as it is in the balloon and satellite measurements |121 |S] where 
energy spectra for different nuclei are obtained up to 10^^ eV). 
Similarly to |411 191j we reconstruct primary energy spectra based on a trans- 
formation of the integral Eq. (jH)) to a parametric equation with unknown spec- 
tral parameters. For that, instead of unknown primary energy spectra of the 
integral Eq. © we consider energy spectra according to the predictions of a 
multi-component model 2111 124j. In this case, unknowns are such model param- 
eters as cut-off and knee energies, spectral indices, scale factors and fractions 
of different components. EAS muon and electron size spectra on the left side 
of the Eq. © are taken from recent KASCADE pubhcations [HI ESI- 
The kernel function of Eq. Q dW/dN^^^ is preliminarily calculated by a Monte- 
Carlo method provided a given (A — A^ir) interaction model of primary nuclei 
and atmosphere. 

Solutions of Eq. Q for unknown spectral parameters at a priori known pri- 
mary spectra are easy to obtain by means of x^-minimization method |9n| I91j. 
It is clear that solutions found that way are partial solutions of the integral 
Eq. Q. However, if the values of the spectral parameters are consistent with 
the predictions of a primary spectrum model calculations [2^, it will confirm 
the use of this model in the investigation energy region. 

We have already done such investigations [91j by testing a multi-component 
model of primary cosmic ray origin with observed EAS size spectra from KAS- 
CADE (1020 g/cm^) and ANI (700 g/cm^) using experiments in 5 different 
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zenith angular intervals. 

Here we consider the inverse problem based on KASCADE EAS electron and 
"truncated" muon size spectra |92| \4:'A\ in 3 zenith angular intervals. The com- 
bined analysis of electron and muon EAS size spectra allows us to determine the 
influence of a nuclear interaction model on the quality of the inverse problem 
solution more accurately. 



9.1 Parametrization of the primary energy spectra 



The energy spectra of primary nuclei according to the multi-component model 
of primary cosmic ray origin j2Ul I24j at energies 10^^ — 10^^ eV are presented 
here in a 3-component form: 



9"^ ( ^ d"^! c d'iS2\ ^EGd'^3 



(4) 



where the /? is a dimensionless normalization parameter, are scale spec- 
tral factors and model parameters 5yi,i=i,2 ai'e the fractions of each component 
{Sa,i + Sa,2 = 1) in a primary flux of nuclei (A). 

The first component (ISM) is derived from explosions of a normal supernova 
into an interstellar medium with expected rigidity-dependent power law spectra 







Ea < EcutZ 
Ea > EcutZ 



(5) 



where the model parameter EcutZ is a cut-off energy of ISM component at Z 
nuclear charge. 

The second component (SW) is a result of explosions of stars into their former 
stellar winds with expected rigidity-dependent power law spectra 



dQ2 
dEA 



Ea'' 

Ek'' 




{Ea/E^ 



-73 



Ea < EkZ 
Ea > EkZ 
Ea > E^JZ 



(6) 



where the model parameter Ef^Z is a knee energy of SW component and Ec,^ ~ 
2.2 • 10^ TeV is a corresponding cut-off energy |2nj . 

The third, extragalactic (EG) component is approximated also by rigidity- 
dependent power law spectra 



^^3 

dEA 



-2.75 

^ank 
p-2.75 

Ea 



E„J,:'° (EA/Eank) 



Ea < EankZ 
Ea > EankZ 



(7) 



The values of model predictions [20] for spectral parameters are: 
71 = 2.75 ± 0.04 , 72 = 2.67 ± 0.03, 73 = 3.07 ± 0.1, Ecut ^ 120 TV, Ek ~ 700 
TV, ^>f£i ~ 0.032 (m^ • sec ■ ster ■ TeV)'^, Eank ^ 6.5 • 10^ TV at factors of 
uncertainty ~ 2. 

In j91j we have already obtained the evaluations of spectral indices of ISM and 
SW components: 71 = 2.78 ib 0.03, 72 = 2.65 it 0.03 as solutions of parametric 
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Eq. © using KASCADE [Hj and |32] EAS size spectra at 5 zenith angular 
intervals. This is in quite good agreement with a similar early test |85| . Since 
these values agreed with the model predictions here we set them fixed. 

9.2 EAS size spectra 

In general, the kernel function of Eq. @ in EAS inverse problems is determined 

by 



dW _ /-^a dG{Eo, A, e) dP{Ne,^) sine 



dNl^ Je, Jo dN,,^ dN^^ Ae 



d9dNe,f, (8) 



where dG/dNe^^ is an expected EAS size electron (truncated muon) spectrum 
at a given observation level for Eq,A,9 parameters of a primary nucleus and 
depends on ^4 — AAir interaction model, A.g = cos^i — cos 6*2, and dP/dN*^^ is 
the error function of measurements |90j . 

We calculated the shower spectra dG/dNe f^ at KASCADE observation level 
(1020 g/cm2) using the CORSIKA6016(NKG) EAS simulation code 03 with 
the QGSJETOl [221 SIBYLL2.1 interaction models. Input parameters 
in simulations were: primary energies Ea = 32, 100, ... 3.2 • 10^ TeV, 4 groups 
of primary nuclei A = 1,4,16,56 and 3 zenith angular intervals similarly to 
KASCADE experimental data j921 143j . Intermediate values are calculated us- 
ing log-linear interpolations. Relative statistical errors of calculated shower 
spectra were less than 3%. 

The expected detectable EAS size spectra dW/dN*^^ were calculated using log- 
Gaussian approximations for dG/dN^^^ and error functions dP/dN*^^ which 
were obtained with high accuracy in the investigated energy range. 
We also investigated the error function of measurements using CORSIKA code 
and for reconstructions of truncated muon sizes we obtained accuracy approx- 
imation: ANfj^/Nfj^ ~ 15/^/Nfj. + 0.05 that exceeds the corresponding values 
of KASCADE approximation 41 . In our calculations we used a more precise 
standard option of CORSIKA6016. 

For the left side of Eq. © we used KASCADE EAS electron and truncated 
muon size spectra from |92| I43j (Fig. 121 symbols). In order to make a compari- 
son and subsequent normalization of our data with direct balloon and satellite 
measurements at about 10^^ energy range we also considered the data from early 
publications of KASCADE |121 with corresponding zenith angular corrections 
(the first 7 symbols for each angular interval in Fig. [2K)- 

9.3 Results 

Minimization of x^(I, P)-functional with a measurement vector 

I^{A//AiV:,,„A//AiV;,^.J (9) 



and a corresponding prediction vector P from the right-hand part of Eq. (j^J 
was carried out at i = 1, ... ,42 energy intervals of EAS electron size spectra, 
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Spectral 
Parameters 



QGSJETOl SIBYLL2.1 3-component 



predictions 



Comments 




2.78±0.03 - 2.75±0.04 

2.65±0.03 - 2.67±0.03 

3.25±0.04 3.25±0.04 3.07±0.1 

200±100 200±100 120-250 

2100±140 1910±150 700-1400 



EH 

3.28±0.07ini] 



0.47±0.04 0.5±0.04 

0.9 0.9 

0.120±0.007 0.106±0.006 

0.089±0.011 0.084±0.010 

0.058±0.007 0.064±0.006 

0.026±0.005 0.035±0.005 



210±60 TV (nU 
1900±100 TV 

0.5-0.8 inu 

0.85-1 HH 
(m^ ■ s ■ sr ■ TeV)~^ 
(m^ ■ s ■ sr ■ TeV)~^ 
(m^ ■ s ■ sr ■ TeV)~^ 
(m^ ■ s ■ sr ■ TeV)~^ 



Table 1: Spectral parameters of reconstructed energy spectra (|4I6() for 4 groups 
of primary nuclei in the framework of a 3-component model of primary cosmic 
ray origin and QGSJET and SIBYLL interaction models based on KASCADE 
date. 



j = 1, . . . , 26 energy intervals of EAS truncated muon size spectra and k = 1,2,3 
zenith angular bins from KASCADE data j92| I43j . However, the combined 
analysis of electron and muon size spectra at x^-minimization requires to include 
in the expected shower spectra 2 additional unknown dimensionless parameters 
rye and rj^ which define a constant bias of each spectrum due to peculiarities of 
interaction model and systematic measurement errors |90j . 
In the first instance, inner spectral parameters (71,2,3, Ecut, Ef^., 6a, ^a) of 
primary spectra (|4I6|) have been evaluated using normalized dimensionless EAS 
size spectra in the measurement vector 



where the scalar parameters Ue = J2i,ki^^ / ^^e,i,k) and Uf^ = J2j,ki^^ / ^^fj.,j,k) 
Corresponding normalizations of spectra were performed for all components of 
the prediction vector Pn at the minimization of the x^(In, Pn)-functional. Af- 
ter evaluation of the inner spectral parameters the values of ijg and r/^ were 
determined using the x^i^^'P) functional by definition Q at fixed inner spec- 
tral parameters. 

Finally, the dimensionless parameter /? in expression (jlj was determined by a 
normalization of the all-particle spectrum obtained ^{dQ/dEA) with JACEE 
data at 100 TeV energy. 

Table 1 contains the values of all spectral parameters which were obtained 
by the method above at QGSJETOl ^ and SIBYLL2.1 [Ml interaction models 
and KASCADE EAS data ^92| I43j . The normalization factors in Table 1 are: 
P = 1.13 ± 0.05 at QGSJET model and /3 = 1.0 ± 0.05 SIBYLL model. 
The energy spectra obtained of different nuclei and the corresponding all- 
particle spectrum in comparison with JACEE 12 and RUNJOB jSj direct 




(10) 
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measurements and EAS data j86i up to energy 3 • 10 eV are presented in 

Fig.m 

The thick hne corresponds to the expected energy spectra of the extra- 
galactic component [2111 ESI EZl EHI • This flux is also included in the all-particle 
spectrum in Fig. ^ 

The expected EAS electron and truncated muon size spectra corresponding to 
our primary spectra obtained with the QGSJET interaction model in compar- 
ison with the KASCADE data [Ml ESI are shown in Fig. ^,h. 

The value of — 1 with 2% uncertainty (the width of a shaded area) of 
expected data (and 4% uncertainty at SIBYLL model). For biases of shower 
spectra (Fig. |2K,b) we obtained the values: r]e = 1.23 it 0.04, rj^ = 1.12 at 
QGSJET model and rje = 1.0 ± 0.01, t]^ = 1.3 ± 0.02 at SIBYLL model. These 
values point out the possible existence of systematic biases both in the EAS 
measurements and in the interaction models. 

In Fig. EJthe expected muon lateral distribution functions (symbols) and corre- 
sponding KASCADE data shaded area) from |^ at different intervals of detected 
truncated muon size are shown. 

In order to test the behavior of primary energy spectra obtained in a larger 
energy range we also calculated the expected hadron energy spectra at mountain 
level for comparison with data of the PION experiment (3200m a.s.l.) jl41 189j 
at 1-7 TeV hadron energies at observation level (700 g/cm^). The calculation 
was carried out using the QGSJET and SIBYLL interaction models and us- 
ing the method introduced above without normalizations (/3 = 1). The best 
agreement of the expected hadron flux at mountain level in the framework of 
the 3-component primary model predictions and PION data was achieved at 
Sa=i,2 = 0.35 ± 0.05 and 5a>i,2 = 0.9 ± 0.1 fractions of component. These 
results are shown in Fig. 0] The effective primary energy range corresponding 
to these data is E'o ~ 10 — 100 TeV/nucleon. 

Thus, predictions of the multi-component model of the cosmic ray origin 
1211 explain the measured KASCADE EAS data in the knee region (£"0 ~ 
10^^ — 10^'' eV) and hadron spectra at mountain level {Eq ~ 10^^ — 10^^ eV) 
with an accuracy of 10-15% in the framework of the QGSJET and SIBYLL 
interaction models respectively. 

The rigidity-dependent behavior of spectra for different primary nuclei is the 
same for the two interaction models. 

The agreement of the expected all particle spectrum and world data in Eq ~ 
10^"^ — 10^^ eV primary energy range displays the presence of the extragalactic 
component of primary cosmic rays in accordance with the 3-component model 
predictions. 



10 The fit and consequences 



Here we have described the actual fit to the data; we did this in two steps, 
first just a direct fit to the KASCADE elemental data as published, as done by 
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2 3 4 5 6 

10 10 10 10 10 

Primary energy (TeV) 

Figure 1: Expected all-particle energy spectrum and energy spectra of different 
nuclei (shaded areas) according to a 3-component model of primary cosmic rays. 
The solid line is the expected (third) extragalactic component and the dashed 
lines correspond to the all-particle spectrum from j21j. The JACEE, RUNJOB, 
CASA, DICE, BLANCA, KASCADE02 data is taken from the review 
KASCADEOl - from 
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Figure 2: KASCADE EAS electron (a) and truncated muon (b) size spectra 
(symbols) at 3 zenith angular bins (from jHS]) and corresponding expected spec- 
tra (shaded areas) in the framework of a 3-component model of primary cosmic 
ray origin and the QGSJET interaction model. 
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Expectation values 
\-g-K --a -2.9-33- - A-4.1.4.5- 

- 3.3-3.7 ■ 4.5-4.9 

t 3.7-4.1 • 4.9-5.3 





40 60 80 100 120 140 160 180 200 

Distance to shower core (m) 



Figure 3: KASCADE muon lateral distribution functions (shaded areas) at 
different truncated muon sizes and the corresponding 3-component model pre- 
dictions with the QGSJET interaction model (symbols). 
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Figure 4: Vertical hadron energy spectrum at mountain level (700 g/cm^) 

obtained by the PION experiment (symbols) and the expected spectra (a shaded 
areas) in the framework of the 3-component primary model and QGSJET and 
SIBYLL interaction models. 
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A. Vasile, and then by an extensive Monte Carlo run to consider all published 
detailed distributions from KASCADE directly, described in the section above 
written by S. Ter-Antonyan. We refer here to the ICRC publications, such as 
[HUEnj. An independent, but similar approach is described by G. Schatz, [5T] : 
recent KASCADE papers are [HllEll^CI- It is gratifying, that a first direct fit, 
as done by A. Vasile, and the extensive fitting done by S. Ter-Antonyan, not 
only agree with each other, but actually agree with the prediction, and with 
the early tests, [55] . 

We summarize, that the KASCADE data as well as a first approximation 
to the HiRes, Akeno and Yakutsk data can be fully fitted by the model with 
parameters very close to what had been predicted. We need to emphasize here, 
that we checked the consistency of the predicted model with the data; we did 
not check the uniqueness of the solution found. We do note that no other model 
has been either inferred nor proposed, which fits all the same data, but it is 
quite conceivable that there is such a model, yet to be worked out or published. 
The approaches by M. Malkov, J. -P. Meyer and L. O'C. Drury, E. Berezho et 
al., and I. Axford et al. all provide a different approach to the quest for the 
origin of cosmic rays. The various approaches may yet converge to a common 
theory. 

Focussing now on the concept originally started some time ago, there are a 
number of important conseqences from the success of the model proposed: 

10.1 A common final state for massive stars 

The are several very interesting consequences from the success of this approach. 
Some of these consequences follow by necessity, which entails if these conse- 
quences can be falsified, then the original model fails. We will identify these 
critical steps below. 

• The abundances in stellar winds near to and beyond 20 solar masses, R,SG 
and WR supergiants dominate the cosmic ray chemical composition - this 
is in contrast to the superbubble concept; however, since the superbubbles 
in fact also have a mixture of wind abundances, these two approaches 
may be consistent in this point with each other. The argument, that the 
abundances of cosmic rays derive from accelerated dust particles is starkly 
different, |65] . 

• The energy in the cosmic ray population of Helium, Carbon, and Oxygen 
is so large as to require a much larger energy content in cosmic rays per 
supernova; note that the RSG and WR stars are the dominant contribu- 
tors, stars which are quite rare. The energy requirement in cosmic rays 
per supernova of this kind is then about 10^^ ergs. 

• If we allow for an overall efficiency of 10 % of putting explosion energy 
into cosmic rays, the explosion energy is then required to be 10^^ ergs for 
such super novae. 

• The most important conclusion from the approach outlined here is that 
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stars near to or above 25 solar masses explode due to the magneto- 
rotational mechanism. We cannot say for sure at this stage at which 
exact mass we have a transition in supernova mechanisms to the magneto- 
rotational mechanism, since the cosmic ray argument is really clear only 
for the WR stars. 

• Knee energy and cutoff energy need to be the same for all supernovae 
contributing. This is easiest if those stars form a "common final state", 
as regards rotation, magnetic fields, mass, explosive energy. 

• There is first evidence for an explosion energy of order 10^^ ergs from 
a very massive star from SN1998bw, [481 169j . That specific supernova 
was also highly asymmetric, with a rotational symmetry quite possible, 
consistent with the mechanism proposed by Bisnovatyi-Kogan. 

10.2 The magneto-rotational mechanism 

The magnetorotational supernova mechanism was suggested by Bisnovatyi- 
Kogan, 29 . The idea was to get the energy for the explosion from the rota- 
tional (gravitational) energy using the magnetic field. Most stars are rotating 
differentially. The magnetic field is "frozen" into the ionized gas of the star. 
Then differential rotation leads to an amplification of the toroidal component of 
the magnetic field. When magnetic pressure becomes comparable with the gas 
pressure a compression wave forms and moves out to the envelope of the star. 
When moving along a steeply decreasing density profile this wave transforms 
to a strong fast MHD shock which pushes part of the envelope of the star to 
the outside, and so forms an explosion. Simulations of the magnetorotational 
mechanism for a magnetized cloud |1U1 166j show that this mechanism is rather 
effective. The simulations of the collapse of a rotating pre-supernova star show 
that after the short collapse stage the star consists of a rapidly rotating core 
and a slowly rotating envelope. At the transition region between the core and 
the envelope the rotation is very differential and in such a situation the in- 
clusion of even a weak initial magnetic field could lead to amplification of the 
toroidal component and so to a magnetorotational explosion. Further work is 
in [SOI El HI]. 

10.3 Imphcations for Gamma Ray Bursts? 

There is an interesting resemblance of this explosion mechanism to the concept 
of a hypernova developed in the context of Gamma Ray Bursts, [ZOIIZS]- Using 
the concept, [7S1 [THI, does indeed suggest an energy scale of also 10^^ erg, 
derived in that specific model as an upper energy limit implied by a fit to 
the fluence distribution, allowing for observational selection effects. The recent 
observations by Schaefer et al., [SO], support the concept that at least some 
GRBs blow up into a stellar wind. 

This leads then one more time to the question of what the relation is between 
Gamma Ray Burst explosions and Supernova explosions, a question explored 
by many. 
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There is much recent work also by C. Wheeler, et al, pi IHl HUl lUllI^ . 

10.4 A bright standard candle for cosmology? 

There is a corollary also for cosmology: If the explosion energy is the same for all 
supernovae above about a zero age main sequence mass of near 25 solar masses, 
as suggested by the arguments above, one may speculate even further: Could 
it be that also the maximum luminosity, or the integral over the luminosity 
time curve is the same for all such supernovae? There is the obvious problem, 
that these explosions are rotationally symmetric only, and so present vastly 
different views to observers at different angles to the symmetry axis. This may 
be resolved by infrared observations, and also by polarization observations, 
just as for active galactic nuclei, for which infrared observations give the best 
approximation for an isotropic emission, and polarization observation do give 
clues for the angle between symmetry axis and line of sight. 

If we could thus derive from observations the maximum luminosity or its 
time integral, we might have a much brighter standard candle for use in cosmol- 
ogy than the supernovae of type la, which derive from the collapse of a white 
dwarf in a stellar binary system, and have a fairly low explosion energy. 

10.5 Transition: Beyond GZK cutoff 

Let us finish with a brief note on the events beyond the Galaxy, the events with 
an energy beyond 310^^ eV. Recent reviews are in, e.g., ^lEl- The latest 
data and their discussions are in PIITI for HiRes, and in [HZl EHl HI IS31 EH EH] 
for AGASA. 

The main dispute at the present time is the apparent discrepancy between 
the AGASA data and the HiRes data, with the HiRes data suggesting a down- 
turn of the overall spectrum in qualitative agreement with the GZK-cutoff con- 
cept, although at a somewhat higher energy, as would be expected from a highly 
inhomogeneous source distribution, such as any subpopulation of galaxies. On 
the other hand, the AGASA data suggest a continuation of the spectrum to 
higher energies, which would require a new component to appear from under 
the lower flux contributors, in perfect agreement with a scenario involving the 
decay of topological defects. The real discrepancy between the two data sets 
is only about 2 sigma, as discussed at the meeting in Paris June 2002. It is 
widely expected that AUGER data will resolve this issue, and then maybe the 
further experiments will be required, such as EUSO and OWL, the European 
and American space missions to detect airshowers at extremely high energies. 

11 Conclusions 

The successful fit of the new KASCADE data, under the constraint to also fit 
the higher energy data from HiRes, AGASA and Yakutsk, with the cosmic ray 
source model originally proposed some time ago, and developed further since 



then, leads to a number of important conclusions. The most important points 
are: 

• The origin of Galactic cosmic rays may be much closer to a full under- 
standing of their origin. Sofar all quantitative tests, which have been 
made, show consistency with the model proposed earlier. But much work 
remains to be done. 

• All stars above a zero age main sequence mass of about 25 solar masses 
converge to a "common final state" . 

• Those stars explode with about 10^^ ergs in a highly anisotropic explo- 
sion, following the mechanism originally proposed by G. Bisnovatyi-Kogan 
more than 30 years ago, which was then based on a broader suggestion 
by Kardashev, involving rotation, magnetic fields and potential energy. 

• These supernovae may constitute a new very bright standard candle, use- 
ful in cosmology, provided we could determine their luminosity integrated 
over 47r from infrared, polarization or other observations. 
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